RIGOROUS WKB FOR FINITE ORDER LINEAR RECURRENCE 
RELATIONS WITH SMOOTH COEFFICIENTS * 



OVIDIU COSTIN AND RODICA COSTINt 
Abstract. We study the e — » behavior of recurrence relations of the type _ %'(fee, t)yk+j = 0, 

t 6 Z (1 fixed). The aj are C°° functions in each variable on / X [0,eo] for a bounded interval / 
and eo > 0. Under certain regularity assumptions we find the asymptotic behavior of the solu- 
tions of such recurrences. In typical cases there exists a fundamental set of solutions in the form 
{exp(e~ 1 i ? m (fce,e))} m=1 j where the functions F m are C°° in each variable on the same domain 
as the aj , showing in particular that the formal perturbation-series solutions are asymptotic to true 
solutions of these recurrences. Some applications are also briefly discussed. 
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1. Introduction. In the present paper we study the asymptotic behavior to all 
orders in e — > of the solutions of one-dimensional recurrence relations of the form 

I 

(1.1) ^2aj(ke,e)y k+j = 

3=0 

which we may interpret as follows: for each fixed k, yk+i is determined from its 
predecessors y k ■ ■ .yt+i-i (this is assumed possible — see condition <|1.6[) below). 

Under some further regularity assumptions we prove that the general solution of 
the recurrence can be piecewise represented as a sum 

i 

(1-2) yk,e = ^2 C I „ iP exp[e~ 1 F m (fce, e)] 

m— 1 

where the functions F m are everywhere smooth with the exception of a small neigh- 
borhood of the points where two characteristic roots (|2.1() cross and where the rep- 
resentation is different (Proposition 2.2) . In particular, a fundamental system of 
solutions can be chosen such that each of them has, for small-e and between cross- 
ings, a WKB-like expansion: 

(1.3) y k - exp (&f>(ke)) (A (ke) + eAi(ke) + ...) 

where <f> is the root of the eikonal equation 

exp(0'(x)) = X(x) 

and A is one of the I roots of the characteristic equation l|2.1|l and the successive 
amplitudes A4 can be determined by perturbation-expansion. For technical reasons, 
we prefer the less familiar notation 11.2(1 . It is essential for our arguments that a 
continuous branch of In Aq can be chosen. 
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One of the applications of the rigorous WKB-approach for discrete schemes is 
in determining the spectrum of large matrices with slowly-varying entries. Such 
problems appear for instance in the continuum limit of the Toda lattice. This system 
, described by the Hamiltonian H = | Ylk=i v\ + Efe=i ex l?{ x k — £fe-i) is completely 
integrable; this can be expressed in terms of the constancy of the spectrum of the 
matrix 



(1.4) 



/ ai h 
bi a 2 b 2 
b 2 a 3 

v 



The spectral problem in case the coefficients are of the form = a(ke), bk = b(ke) 
with a and b smooth and satisfying some regularity conditions leads to a recurrence 
of the type that is solved asymptotically by the methods described below [1]. 



In the number I of steps of the recurrence is fixed, e > is a small parameter 
and k G 7L is such that ke € I where 7 c R is a compact interval. Some initial or 
boundary conditions are assumed. 

The coefficients ao(x, e) . . . a/(x, e) : 7 x [0, e ] — >ff are assumed C°° in x and in e 
in some domain 7 x [0, eo]. We also require the existence of a uniformly asymptotic 
series for a,j\ for any t G TV, 



(1.5) \a 3 {x,e)~Y, a ^ x y \< M J^ +1 

where the functions aj iS are C°° in x (for instance aj G C°°(7 x [0, eo]). 
We are also imposing the nonsingularity condition 

(1.6) inf{|a (.T,0)|} > and inf {\ai{x, 0)|} > 

We begin by giving some simple examples and deriving heuristically their small- 
e behavior. The contents of the paper will subsequently make the given solutions 
rigorous. 

a. Consider the one-step recurrence relation 

(1.7) Vk+ i = e /(fee V 

where y = 1, < k < e _1 and / is a C°° function on [0,1]. It has the explicit 
solution 

fc-i 

(1-8) 2/fc=e X p((£/(je)) 

3=0 

When / = /o is constant, = exp(e~ 1 /ofce) and this should also be the order 
of magnitude of yu for a general smooth f when e is very small; we then try a formal 
asymptotic solution 

(1.9) y fe ~ exp(e- 1 $ (fce) + $i(fce) + e $ 2 (/ce) + ...) 
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in l |1.7f l . To be consistent with (jl.7(l we must have 

oo oo 

(1.10) exp[^ e m -^ m {ke + e)] ~ cxp[/(fce) + £ e^^ke)], e -> 

m— m— 

Expanding the exponent in the LHS of (| 1 - 1 U|> in a Taylor series around ke and then 
identifying the corresponding powers of e in ll.lOfl we obtain 

(1.11) $0 = /, ^ = ~\f ••• &i = C i /°' ) .-- 

where the constants Cj are seen to be /-independent and can thus be determined by 
choosing a particular / for which the sum in l|1.8[) can be explicitly evaluated. Taking, 
e.g., f(x) = exp(x) this sum is 

oo 

(1.12) (exp(x) - l)(exp(e) - l)" 1 ~ e'^exp^) - 1) £ j^r^ , e -> 

where B n are the Bernoulli numbers. It follows that Cj = Bj (Baj+i = 0) and thus, 
for any smooth / 

(1-13) Hm--X^ n f (n) (W, e^O 

j=Q n=0 ^ >' J ° 

which is, of course, the Euler-Maclaurin summation formula. Proposition 12 . II below 
justifies our derivation of Q1.13|l . 

b. Analogously, we can obtain easily the asymptotic behavior of special functions 
from their generating recurrence relation. For instance, the recurrence relation 

(1.14) y k+1 +y k -i = 2{l + ke)y k 

has, for a fundamental set of solutions the Bessel functions Jk+e- 1 ( e_1 ) an d Y k+li -i (e _1 ).| 

Let x — ke, v = k + e _1 , a = e~ x v~ x . We fix p < | and for for x > e p we take 
y of the form l|1.9fl determining the successive terms by substituting the formal series 
in l|1.14|l . To leading order 

(1.15) e *iW +e -*;W =2 (i +I ) 
and we can choose two independent solutions 

(1.16) e _1 $o ; ± = ±v(a - tanh(a)) (a = (cosh) _1 (l/a)) 
For the next order we obtain: 

(1.17) $' 1 . ± /$ 1;± = -i^'. ± coth($') 

which can, again be integrated explicitly to give exp($ 1; -|-) = const sinh(a) -1 / 2 . That 
is, the asymptotic behavior is 

(l.l^tanh(a)) _1/2 (Aexp[i/(a - tanh(a)) + ...]+ B exp[— v{a — tanh(a)) + ...]) 
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and we get the familiar expressions in the theory of Bessel functions (see [2], [3]). All 
the successive orders are obtained easily in the same way. For x < — e p we obtain 
from the same equations (|1.15|l . (|1.17[) similar expressions but with trigonometric 
instead of hyperbolic functions. For small values of \x\, \x\ H e 2 / 3 the above 
asymptotic series becomes singular. The appropriate WKB-like series in this region is 
in powers of e 1 / 3 and the coefficients will be smooth functions of e 1 / 3 ^, but otherwise 
the calculation can be done explicitly in the same way; there is a region of overlap 
where both asymptotic series are valid, namely e 2 / 3 < |x| < e 1 / 2 and where the series 
can be matched. 

Proposition 12 . 21 can be used to make the above approach rigorous. 

c. Finally, let q be a smooth function and consider the Cauchy problem 

(1.19) y" = g(x)y, 2/(0) = 0, y'(0) = l 

Assume for simplicity that q : [0, 1] — * 1R + and consider the associated Euler scheme 
(1-20) Vk+e-.e + Vk-e-,e = (2 + e 2 q(ke))y kt e, y . e = 0, y 1;e = e 

To characterize yk-e for small e we substitute 

oo 

(1.21) y k ;e ~ ^f m S m (fc) 

m— 

in (|1.2U|) . Note that if So ^ we obtain this type of series from l|1.9|l . when $o = 0, 
by expanding the exponential. 

The substitution leads to the equations: 

(1-22) 5" l = q(x)3 m -2 £ 

l<s;2s<m ^ >' 

with the initial conditions 

m CT (fc+l)/Q\ 

(1-23) 3 m (0) = 0; 3^(0) = " E ^TT)T 

In particular, it follows that the scheme converges to the solution of the given Cauchy 
problem (as it should) and because, as it is easily seen, Si = the error |yjc+ e — S (fce)| 
is 0(e 2 )- Proposition 12.31 applies to this example. 

2. Main results. In this section we consider the problem (| 1 . 1 1) and the assump- 
tions following it and give conditions under which the exact solutions of the recurrence 
have asymptotic series to all orders in e. 

Proposition 2.1. Let Ai(a;), . . . A; (a;) be the roots of the characteristic polynomial 

i 

(2.1) ^a J (x,0)A^=0 

and assume that they are simple throughout I: 

(2.2) inf{| X m {x)-X n {x) \}>0 ifm^n. 

xEI 



(consequently, we will choose X m to be C°° functions). Suppose also that the interval 
I is a finite union of (closed) intervals such that in each one of them the ordering of 
the moduli of the roots does not change, i.e. 



I = I p and Vp < P 3(ii . . .ii) such that 
P =i 

(2.3) \X h (x)\ < \X i2 (x)\ < ... < \X k (x)\ for all x e I p . 

where (i\ . . .ii) is a permutation of (1 . . .1). Then there exists e < eo and I func- 
tions F%(x, e) . . . Fi(x, e) on I x[0, e'], C°° in each variable, such that for each e, 
{exp(e _1 i 7 ' m (fce, e))} m= i...i form a fundamental set of solutions of the recurrence re- 
lation in the sense that for any solution y^ e of there exist constants 

C| , . . . Cj 1 ^ such that in each L t 



(2-4) y k , e = C ™ exp[e- x F m (fce, e)] 

771 = 1 

Remark 1. In particular, this means that for small e 

oo 

(2.5) F m (a;,e)^^$ m , s (x)e s 

s=0 

where the $> m ,s o,re smooth; they can be computed from by usual perturbation 

expansions in e. For example the first term is gotten from the eikonal equation, 

(2-6) expffc^x)) = X m {x) 

giving the connection between the functions F m and the roots X m of the characteristic 
polynomial. The second term is obtained from 

(or , , > _ Ej =0 [i 2 ^n,o(x)72^,o(^) + exp(<fy {x))Aj, iQc)] 

E J= o^j,o(^)exp($;- (x)) 

and so on. 

Note also that if not all the functions F m have the same real part, then there exist 
solutions of the recurrence which are exponentially small relative to the dominant 
ones and which will therefore be unstable in the sense that a small "generic" pertur- 
bation of the initial condition will change completely the behavior of these solutions; 
a perturbation series not involving terms beyond all orders will only see the compo- 
nent of the solution along the dominant directions. However the relative size of the 
solutions could change with x and then all the coefficients in the expansion 12.4fl could 
be important for matching different regions. 

We now address the question of the asymptotic behavior of the solutions of the re- 
currence when two characteristic roots cross. The previous Proposition is generalized 
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below to the case when two of the roots of the characteristic polynomial become equal 
at some point in / provided the roots do not coalesce too quickly (condition J2]SJ). 
In this case a fundamental set of solutions has a more complicated structure. Not too 
close to the crossing point, a solution is still a linear combination of the form i|2.4|) but 
the coefficients Cm) can change at the crossing (Stokes phenomenon). Very close of 
the crossing, the coalescing roots bring in the asymptotic expression of the solution, 
series in noninteger powers of e. This is in some sense the discrete counterpart of the 
turning point-behavior of the solution of a differential equation depending on a small 
parameter. 

Consider a subinterval Io of I such that two roots of the characteristic polynomial 
cross once in Iq (say, at x = 0) and except for this, the ordering of the moduli of the 
characteristic roots is constant in Iq. The crossing is assumed to be generic: 

(2.8) |A m (x) — A n (a;)| > const ( for to 7^ n) 

To avoid excessive branching of the discussion and formuli, we also assume that 
the coalescing roots are complex- conjugate for negative x and real-valued for x posi- 
tive. The general case is treated in a very similar way. 

Fix two constants | < f3 < a < |. Then, 

Proposition 2.2. i) For \k\ > e^ 13 the general solution of the recurrence can be 
written in the form 

1 

(2.9) 2 / fc = ^C,cxp(e- 1 F,(fce,e)) 

i=i 

with Fj as in Lemma \3.°2\ The constants Cj depend, in general on the sign of k. 

ii) For \k\ < e~ a a fundamental set of solutions can be chosen such that: 
I — 2 solutions are of the form exp(e _1 i^ (fee, e)) and two special solutions of the form 

exp(F±(fce 1 / 3 ,e 1 / 3 )) 

where the functions Fj and F± are smooth and exp(F±(x, 0)) —Ai(Q x)±Bi(<3 x) 
where Ai, Bi are the Airy functions (for the value of see equation \4-14[ ) ■ 

Moreover there is a particular solution of the recurrence which has the behavior 



y k ~ Ai(G fce 1/3 )(l + e 1 / 3 J 4 1 (6 fce 1/3 ) + ..) 

for large k < e~ a 

The representations in (i) and (ii) are simultaneously valid in the region e^ 13 < 
\k\ < e~ a where the asymptotic series can be matched. 

Note. A similar result can be proven if condition 1)2.8(1 is replaced by 

(2.10) \\ m {x) - \ n (x)\ > const x const ' 

for some const > 0. Another case which is interesting for schemes that converge to 
differential equations for small e is covered by the following proposition, in which we 
assume that the roots of the complete characteristic polynomial: 

/ 

(2.11) '^a J (x,e)X j , m=l...l 

3=0 
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are nondegenerate in a higher order in e. 

PROPOSITION 2.3. Assume that the roots Xi(x, e) . . . Xi(x, e) of \2.11\) satisfy the 
estimates 

(2.12) \X m (x,e)-l\=e q (Q m (x)+o(e)), m<l 
for some geff where the smooth functions Qi verify 

inf Q m (x) > 0, inf \Q m (x) - Q n (x)\ > (m ^ n) m,n < I 

Then, the conclusion of Provosition WA\ holds and, furthermore , for any formal series 
solution of the recurrence relation there exists a true solution which is asymptotic 

to it. 

In this particular case, since, as we shall see, F m (x,0) = 0, it is more natural to 
represent the formal solutions as power series: 

oo 

(2.13) r = Y / * l (xV 

where the 4^ are smooth and subject to the condition 

i 

(2.14) ^a i (a; ! e)r s (x + je;e) = o{e s ) Vs, x E I. 
j=o 

where T s (x; e) = Yli=o ^i( x ) el - The series (I2.13f) and those appearing in the previous 
WKB expansions are usually divergent and one could imagine that by iterating the 
recurrence the small error appearing in the local condition (I2.14|) could quickly reach 
0(1). Under the given restrictions, however, Proposition 12.31 guarantees that T s is 
indeed an o(e s ) approximation to a genuine solution. 

It can be now checked without difficulty that Propositions 12 . II 12 . 21 and 12 . 31 apply 
to the examples a), b) and c), respectively. 

The layout of the paper is as follows: in Section 3 we prove our main results and 
in Section 4 we discuss some further applications of these results. 

3. Proof of the results. To prove Proposition 12 . II we show (Lemma I3.1fl the 
existence of I formal series solutions of the form (|1.9|) to the recurrence. 

We then show (Xemma |3.2[) that the the proposition is true if P = 1 (cf.J25J)). The 
proof is by induction on the order I of the recurrence. First, we choose the particular 
formal solution corresponding (cf. I|3.2|l ) to the root with maximum modulus - 
which gives the "stable" direction — and show that there is a true solution with this 
asymptotic behavior. Next we use this particular solution to decrease the order of the 
recurrence by 1. 

3.1 The functions $ mjS of the following lemma will turn out to be the functions 
giving the asymptotic expansions H2.5JI and can be obtained by requiring that 1)1. ljl is 
satisfied in all orders in e by the formal solution y m — exp(e _1 J2 S ^m,s{x)e s ). 

Lemma 3.1. For each m = 1 . . . I there exists a sequence {^ m ,s} s& ]]\[ of functions 
in C°°{I) such that 

l s 

(3.1) exp(- e - 1 $ m , (fce))^ a ,(fce, e )exp(e- 1 ^ e'$ m , t ((fc + j)e))) = 0{e s+1 ) 

3=0 t=0 
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and 



(3.2) exp« i0 (z)) = X m (x) 

for e < eg, s G IN and ke G / 

The proof of Lemma 1 is by induction on the expansion order s. In view of (|1.5f) 
and 1(1.6(1 we can define for each characteristic root X m {x) (i.e. root of eq. 1(2. lp l. a 
function $ m , G C°°(I) such that (|3~2j) holds. 

It is then straightforward to show that exp(e _1 $ m .o) verifies (11.1(1 to 0(e) exp(e~ 1 $ m .o(a ; ))| 
so that (|3.1() holds for s = 0. Assuming now that Q m .o, $ m ,i ■ • ■ ®m,s are already 
defined so that for all s < so 1)3.1(1 is verified, one can easily check that for any 
* G C°°(I), 

exp[-e" 1 $ m:0 (/ce)] 
J2 expj^ 1 *"»,*((* + fry + e S0+1 *((^ + J»K(^ e) - 

j=0 i=0 

i 

(3.3) e*° +1 ]T j aq){i^ i0 (fce)}o i (fce ) 0) + H so (ke, e)] + 0(e s °+ 2 ) 

3=0 

where H sa is a smooth function. 
Since by {HI and l|23| 

i 

(3.4) inf |^ja i (x,0)A4(x)| >0 

one can define a smooth function ^(x) = (fr So +i(x) such that the term in square 
brackets in the RHS of l|3.3|) vanishes. Oli 

We note at this point that the series ^m,s(x)e s is, usually, not convergent 

and there does not yet follow the existence of a solution asymptotic to it. 

Now we address the question of existence of true solutions of the recurrence hav- 
ing the precribed asymptotic behavior. In what follows, we shall understand by a 
formal solution an expression Y = cxp(e _1 ^m,s(x))e s satisfying the conclu- 

sions of Lemma 13. II Given $> m ,o, the & m ,s are uniquely determined up to integration 
constants. 

Assume first P — 1 (cf. 1(2.3(1 ). Relabelling if necessary, we assume that for 
m<n, |A m (x)| < |A„(a;_)| on /. 

Lemma 3.2. Let Y be a formal solution of (m G {1 . . .1} fixed). There 

exists a sequence {Y mt k-e}k,e such that for any e < eo, Y m ,k-e is a solution of 
for ke G / having T m as an asymptotic series in the sense that there is a sequence of 
positive constants {C s } s such that 

s 

(3.5) |F m , fe;e exp[- e - 1 ^$ mit (fce)e t ] - 1| < C s e s+1 for all s G W 

t=o 



s 



Remark 2. The previous lemma can be restated as follows: for each m=l...l 
there is a function F m (x,e) : / x [0, eo] — > (E , smooth in each variable, such that 
exp(e F m (ke, e)) is a solution of the recurrence $1.1}) for every e, ke E I and, as 
e -► 0, 

OO 

(3.6) F m (a;,e)~^$ m , 4 (x)e* 

t=o 

This remark follows easily from a classical result ( see e.g. [4], page 33 and [5]) stating 
that for any sequence of numbers there is a smooth function having that sequence for 
its derivatives at the origin and from the easily proven fact that for each sequence 
{(x n , <Xn)} n c]N > x n \ and a n n k — > Vfc one can construct a C°° "interpolation" 
function f such that f( k '(0) — Vfc and f(x n ) = a n for all n. 

Some of the estimates that we need for proving Lemma 2 become to a certain 
extent easier by the remark that a global shift in the estimates defining an asymptotic 
series is unimportant: 

Remark 3. 

Let F be a function such that for a fixed sq > and any s in IN 

S+SQ 

(3.7) limsupe- s - 1 |F(e) - V C t e*| = D s 
Then 

s 

(3.8) limsup e - s - 1 |F(e)-V C t e*| <D S + \C S+1 \ 

Thus, to prove (|3.5() we need only show that for some fixed sq, (we drop the 
subscript e to ease the notation) 

s 

(3.9) \Y m . k expt-e- 1 ^ $ m , t (fce)e*) - 1| < C s e s+1 - Sa 

t=o 

Remark 4. In view of condition it is easy to check that, for small e, the 

solutions exp(e _1 i r ' m (fce, e)) are linearly independent, thus forming a fundamental set 
of solutions on I . 

Proof of Lemma 13 . 21 The proof is by induction on the order I of the recurrence. 

Step 1. For any I, in the same conditions as in Lemma \'3.1\ the estimate $3.5\) 
holds for m=l (recall that \\{x) is the largest in absolute value). First we show that 
if a solution is asymptotic to the formal solution corresponding to Ai at the left end of 
/ it remains asymptotic to it throughout /. We can assume without loss of generality 
that the left end of I is at x = 0; at this point we choose appropriate initial conditions: 

let T] p , p = 1 . . . I be any functions such that as e — * 

oo 

(3.10) r] p (e) ~ exp(e- 1 £ $ M (pe)e') 



t=o 



(see Remark 2) and let 



(3.11) Y £ , s (x) = expie- 1 V $i,t(x)e*) 



t=o 
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Let also Y\ k be the solution of satisfying the initial conditions 

Yi tP = T]p(e) , p = 1...L It is natural to rescale the recurrence relative to its 
approximate solution: let 

(3.12) C k;s =Y ltk /Y e . iS (ke) 
Then the recurrence relation for C k]s can be written 

l 



(3.13) 



^2a jtS (ke, e)C k +j- s = 
j=o 



where aj ;S (fce,e) = cij(ke, e)Y e - s ((k + j)e)/Y e - s (ke). It can be seen that (|3.13() is of 
the same type as ljl.lt and that, for e small enough, it satisfies the corresponding 
assumptions (|1.6(l . 12.2(1 and 1(2. 3[) on /. 

Also, max{|A m (a;)|; x E I, m = 1 . . . 1} = 1. Now, 1(3. 9|) means that for some fixed 

so, 

(3.14) \C k;s - 1| < const s e s - Sa+1 

From the definition of Y t . s it follows that 



(3.15) 



5>, s (fe,e) = 0(e s ) 



Rewriting the recursion relation (|3.13|l in matrix form Cfc + i = A k C kl where 

/ -a;-i(fce,i;) -q;- 2 (fce,e) -ao(fce,e) \ 

1 ... 

(3.16) 



V 



and rewriting also ((3.15(1 as 
(3.17) 







A k l = 1 + e s E fe 





o / 



where 1, = 1 and ||Efc|| < const uniformly in fc, e we have 



(3.18) 



k 

C k = l + e s ^2 A k A k -x . . . A j+ iE k 

3 = 1 



Step 1 is completed by showing (|3.14l) (and thus 1(3.9(1 ). which follows from the stability 
lemma below. 

Lemma 3.3. Let A k be a family of matrices of the form fcl.lb]) where hj : I x 
[0, eo] — >(F (I is an interval) as in Lemma \ti.'i\ 

Suppose further that the roots A m (fce,e) of the polynomial 



(3.19) 



5^5^,6)^ = 0, 

.7=0 
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m = 1 . . . I 



satisfy 

(3.20) sup{|A m (fce, e)|} < 1 + const e, m=l..l 

keel 

and that the condition corresponding to \2. ZH) is fulfilled on I. Then there is an e- 
independent constant C such that 



(3.21) \\A k A k ^...A ]+1 \\ < l + C\k-j\e 

Proof of Lemma 3. The eigenvalues of the matrix H3.16(l are the A m (fce, e) and 
the corresponding eigenvectors matrix is (Ak)ij = Xj(ke,e) l ~ 1 . We then write the 
product on the LHS of (|3.21|) as 



(3.22) A k D k A^A k -iD k ^ . . . A ]+1 D ]+1 Aj^ 
where 

(3.23) D p = diag({\ m (pe, e)} TO= i...;) 
and the proof follows from (|3.2()|l and the estimate 

(3.24) HAj^Ap-xH < 1 + const e 

which can be checked, for instance, using the following explicit formula, whose ele- 
mentary proof we omit: 

Remark 5. Let X and Y be two nonsingular Vandermonde-type matrices Xi_j = 
x„~\ = y l r\ i,j = l...l. Then, 

— l\r\ _ TT Vj 



(3.25) (X-'Ykj = ]J 



Step 2 . The conclusion of Lemma 13.21 for I = 1 follows from step 1 . 

Now we assume that the conclusion of the lemma holds for all recurrences of order 
less than I and prove it for order /, by reduction to the I — 1 case. In view of the 
first step, we know that to |Ai| there corresponds a true solution Y\ for which the 
asymptotic behavior is the formal solution Y\. We shall use this solution to reduce 
the order of the recurrence by one. Let 

(3.26) C k = y k /Y hk 

The recurrence relation for C k is then of the form l|3.13l) where now 

(3.27) a,j(ke,e) = aj{ke, e)Y 1:k +j/Y 1:k 

and obviously, the asymptotic behavior to all orders is the same as if we had made the 
rescaling with respect to a formal solution. The point is that now, instead of l|3.15|l 
we have 

i 

(3.28) ^a i (fce,e)=0 

j=o 

li 



so that the y = 1 is an actual solution. To use this fact, let dj. = Ck+i — Ck- We get, 



i-i 

(3.29) y £ / b s {ke,e)d k+s = 

s=0 

where b s (ke,e) = I^- =s+1 Oj(fce, e). 

The characteristic equation for l|3.29|l can be written as 

/ j-l 
^a,(x,0)^A s =0 

j=l s=0 

or, for A ^ 1, as it easily follows from (|3.28|) . 

I 

(3.30) ^a i (x,0)A J =0 

We first check that the new recurrence satisfies the hypothesis of the Lemma. But 
this is easy since the new coefficients are finite combinations of a,j and in particular 

&o = Sj=i % = — °o an< i l = a d an( i m v i ew °f H3.30fl the same arguments as in 
step one apply to see that the characteristic roots have the required properties. Then 
we want to check that we have the required number of appropriate formal solutions. 
This is also straightforward because we can derive them from the formal solutions of 
the original equation. Indeed, 

(3.31) d k = Vk+i/Yi,k+i - Vk/Y\,k 

If we substitute for y a formal solution Y m we obtain the formal expression : 

oc 

d k = expK 1 X>m,t((fc + l)e) ~ $!,*((* + l)e)]e*} 

4=0 

oo 

(3.32) - expje" 1 ^[$ m:t (fce)) - $ 1)t (fce)]e*} = 

t=o 

exp {e _1 ($ m (x) — $i(a;)) + series} [A m (x)/Ai(a;)](l + e x series] 

which, since A m (x)/Ai(:r) does not vanish can be written as an exponential of a 
formal series, the form required by our arguments: 

oc 

(3.33) exp( e - 1 ^A m:t (fc e )) 

t=o 

and then the expressions (|3.33l) for m = 2 . . . I are formal solutions for the recurrence 
of order I — 1 i|3.29|l — because, by construction, l|3.32|l are; it follows, by the induction 
hypothesis that there exist true solutions of l|3.29|l of the form 

(3.34) d m .u = exp(e _1 D m (fee, e)) 
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where -D m (-, •) are smooth functions having the asymptotic behavior given by (|3.33|) . 



Step 3. To complete the proof of Lemma 2 it remains only to check that 

k 

(3.35) exp (e^F^ke, e)) ^ exp(e~ 1 D m (pe, e)) 

has the asymptotic behavior needed for the original recurrence, i.e. 



(3.36) 



exp(e- 1 2^($„ M (fce)e*)) 



We let k\ (^2) be the left (right, respectively) end of the interval. Both k\ and ki 
might depend on e. By the definition of the dk- m we have 



With the choice 



i—ki 



we get the particular solution 



k 2 



from which, referring to the definition of the C m -k, we get a solution of the I — 1 
recurrence in the form (the choice of the sign will become clear later) 



fe 2 



i=k 



whose the asymptotic behavior is given by the formal solution Y m . Indeed, for any 
fixed large s we have, 



Ym:k — —Yl:k 



£exp £ -^A f (^ (l + O^ 1 )) 



,i—k 



-Yi- 



^ exp e 1 + e)e* - $ 1;t (ie + e)e* 



1.3 



exp e- 1 J2 * TOi t(ie)e* - $i;t(ie)e' (l + O^ 1 )) 



t=o 



exple 1 ^$ m;t (fce)e* - $i ;t (fce)e* -exp e 1 ^ $ m; t(fc 2 ) - $i ; t(fe)e t 



t=o 



t=o 



(3.37) 



(&2 — k) max 



exp e > -A/(/r )' 

\ t=0 



0(e s 



Because, by assumption, Ai has the largest modulus, 5i($ m; o(fce) — $i ; o(fce)) is 
nonincreasing in k in the given region. Therefore l|3.37|l equals, 



Y, 



1 + const max { | exp ($m-x(ke) 
k 



^ ul (ke))\}o(e s - 2 ) ) = 



(3.38) 



Y m , k (1 + o(e s - 2 )) 



UL2 



The proof of Proposition ^ . 21 follows essentially the same steps but is, as expected, 
more involved in the regions of near-breakdown of the asymptotic series. The details 
are given in the next section. 

The proof of Proposition ^. 31 is very easy, using an estimate of the form l|3.21|) for 
the matrices corresponding to the original recurrence, estimate which is straightfor- 
ward to obtain from Remark 5 and the hypothesis of the proposition. 

4. Proof of Proposition 2.2. We assume at first that the crossing occurs be- 
tween the two largest characteristic roots and explain at the end of the proof how the 
general case is reduced to this one. 

The layout is as follows. We first study the small region around the crossing 
point (the interior region) where exp(e _1 ^i,2;t(x)e t ) fail to be formal solutions 
(and the series occuring at the exponent cease to be asymptotic series). The new 
formal solutions are to leading order combinations of Airy functions. Their formal 
properties (domain of asymptoticity, growth in x) are examined. Next we show that 
there exist true solutions of the recurrence that are asymptotic to them. It is also 
shown that there exists a particular true solution which is asymptotic, to leading 
order, to the function Ai(cce -2 / 3 ) and which is important for the matching problem 
(it gives the exponentially decaying formal solution). 

We then show that the formal solutions coming from the exterior region continue 
to represent correctly the solutions of the recurrence far enough into the interior region 
(down to |x| ~ e 2 / 3 ) to allow for matching with the interior ones, which are valid up 
toM^e 1 ^. 

A. The interior region of the crossing interval 

This is the region \x\ <C ; for definiteness we fix a G (|, |) and take it to be 



D a = {k : |*| < e- a } 
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or, in terms of £ := fee 1 / 3 := kS which is, for reasons that will become clear later, the 
natural variable in this region, 



(4.1) D a = {£ : |£| < 6 1 - 3 *} 

The basic steps of the proof of existence of solutions with given asymptotics are 
the same as for Lemma 13.21 We will first obtain a solution corresponding to (one 
of the two) largest eigenvalues and with it reduce the problem to a lower order, 
nondegenerate recurrence. 

Because in the variable £ D a is unbounded, a slight extension of Lemma [3.21 is 
needed for the interior region. We now allow the interval / in Lemma l3~2~l to be of the 
form 14.111 but strengthen the other hypothesis. In order to make the correspondence 
with Lemma [3.21 note that £ plays the role of x and <5 is the counterpart of e. We 
require the same conditions as in Lemma 13.21 and in addition, 

a) Oj(£, 6) are assumed to have asymptotic series J2t a j,t(0^ t valid throughout 
the region D a which are smooth in the sense that the all their formal derivatives with 

respect to £, J2t a j!t\0^, exist and 1^(01 < const 5 const ' 1 

b) The roots of the characteristic polynomial 

I 

(4.2) ^ aj (£,0)A(^"=0 
are nondegenerate: 

(4.3) infD a \\ m {0 - A„(£)l > const > (m ^ n) 
and the polynomial itself is nondegenerate in the sense 

(4.4) infz,^ I \a Q (k8)\, \ai(kS)\, -, — ^-^ \ > const > 



c) 



\X m ((k + l)S,5)-\ m (k5,S)\ < 

ft + 1 



in D a where A m ((fc + 1)5, S) are the roots of the complete polynomial J2j=o a j(£,i ^)^' 
Lemma 4.1. Under these assumptions for any formal solution of \1. 1)) of the 
form 

oo 

m=0 

where the exponent is assumed to be a smooth asymptotic series in the sense defined 
in a), there exists a true solution of the recurrence which is asymptotic to it in D a . 

Proof. The proof follows closely the proof of Lemma l3~2l We only emphasize the 
differences: For the recurrence (|3.13|) we have also to verify condition c) : 



X m ((k + l)S, 6)-~X m (kS, 6) 



X m ((k + 1)6,6) X m (k8,6) 



Xi({k + 1)8,8) Ai(JW,<5) 



< 



15 



const (|fc| + 1) 1 

The equivalent of Lemma 13.31 states now that there is an <5-independent constant C 
such that 



(4.5) ||A fc A fc _ 1 ...A 3+1 || < (l + C|fc-i| const ) 

Indeed, by the Remark 5 a diagonal term of the matrix T := A^ 1 A p _i is of the 
form 



rp tt X m ((p-l)S,S)~X n (p6,6) _ 

i m , m - Li x m (ps,s)-x n (ps,s) ~ 

««) =n0 + ^r^ , ' ) )° 1+o '" t|+1) -» 

by c). Similarly, the moduli of the nondiagonal terms are seen to be less than 
const/(|fc| + 1). Therefore T = I + R where the ||i?|| is 0((\k\ + l)" 1 ), hence the 
inequality H4.5(l follows. 

The last part of the proof of Lemma 13.31 applies here without any significant 
change. Dl4.i 

The reduction to the nondegenerate case. 

We consider the initial recurrence in the neighborhood of a crossing point, say x=0 
where Ai(0) = A2(0) = 1 (the value at zero can be chosen through a trivial global 
rescaling of the recurrence). It is convenient to consider rescaled variables S — e 1 / 3 
and £ = kS. In these variables, the coefficients aj(x,e) — cij(£5 2 ,6 3 ) have smooth 
asymptotic series in 6 in D a which are in fact obtained through series expansion in x 
from (|1.5ll : 



(4.7) u,^.^)-y^l>,.A^ 

s>0 

where 

P jfi = a,(0,0); P^O = 0; P ja fc) = (D x a 3 )(0,0)^; 

and in general P, jS (£) are polynomials in £ of degree at most s /2 for s even and (s— 3) /2 
if s is odd. To avoid complicating the notation we write aj(^,S) = aj(£,5 2 ,S 3 ). We 
have first to find formal solutions for this new recurrence. 

Lemma 4.2. There exist I linearly independent formal solutions in D a , of the 
form 

(4.8) expJr^-iVttO^ 

where ^ m ,t{C) are smooth in £ and satisfy the estimates 

16 



\^m,t(0\ < constmj + consif m t \^\ 



4+1 



This means in particular that the domain of formal validity of the power series 
is then £-/ 2 8 <C 1 i.e., x <C 1. The domain in which it is actually asymptotic to the 
solution is however much smaller (x <C y/e) as we shall see. 

Proof of Lemma 4.2. 

The formal solutions corresponding to the nondegenerate roots give rise automat- 
ically to acceptable formal solutions in the new variables £, 8. Indeed, 

(4.9) expfe- 1 ^^^^ = 



exp u- 3 $„ M (o)+5- 1 * { &m a s 



^( s ) (r\\c s ^2s+3t— 2 

s,i>0 



The term in 8 3 is merely a multiplicative constant so it can be dropped and we are 
left with a formal solution of the form 



(4-10) exp\^- 1 J2xm,t(£)S t j 

where the Xm,t(Q are in fact polynomials in £ of degree < t/2 + 1. 

For m = 1, 2 it is more convenient to write first the possible formal series solutions 
for the equation, in the form: 



(4.11) X>(0* 



t=o 



and then show that we can write them in the form l|4.10|l . 
Substituting (|4.11|l in the recurrence we get 

I oc 
j=0 t=0 

The term of order s in (|4.11|) is gotten by differentiating the auxilliary equation 



i 

i>(e,*)x(f+jA*) = o 
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s times with respect to 5. s times with respect to 5. We get (see (|4.7H ) 



i 



EE p -*(£) (^+^)*x 

j=0 t=0 







delta— 

which after expansion, change of order of summation and use of (|4.7[) gives, 

s— 3 s 



I 

E 

j=0 



f s \ E E 



t 

tj \<7 



(4.12) iW0 J Q)ex.-2 + a J -(0 > 0)j 2 x.-2] =0 

It follows that xo is obtained as a solution of the homogeneous Airy equation 

(4.13) x"(0 = e 3 £x(0 

where 



(4.14) 



£,=0^(0,0) 



and that, given Xo---Xs-3i we get Xs-2 as a solution of an inhomogeneous Airy equation 
of the form 



x"(0 = e 3 £x(0 + #(0 

where i?(£) is a linear combination of higher derivatives of xo.-Xs-3- To avoid cumber- 
some notations, we shall assume in the following that 8 is one. We can check that the 
assumption of genericity (| Ai(ar) — A2(x)| > const y/x) implies E J= o D x a,j(0, 0) ^ 0. 
We shall assume for dehniteness that it is negative. 

It follows by an obvious induction that the \s are smooth. Now we show that 
they satisfy the inequalities stated in the Lemma 4.2. 

Remark Consider the inhomogeneous Airy equation /"(£) = £/(£) + R(0 and 
assume 

R(0 ~ £ p exp(2,4/3£ 3 / 2 ) with A = ±1 for x — > oo and A = ±i at — oo . Then 
/(C) ~ exp(2yl/3£ 3 / 2 ). This estimate follows immediately from the explicit form 
of the solution: 



/(0=Ai($)/ R(t)Bi(t) dt - Bi(g) J R(t)Ai(t)dt 

At this point we can show by induction that the solution x«(£) grows at most 
like exp(2A/3C 3/2 )r /2 - 

So we assume that this holds for s < n and we show that it 
is true for n+1. Using the remark and l|4.12[) the induction step is: with p n = n/2 



max 

0<er<n;CT<t<n+3 



i(n+l-t) 



< Pn+1 + 1 



which is straightforward. 
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Finally, we argue that there are two linear independent formal solutions of this 
type that can be written in the form (|4.1U|I which is convenient for our approach. 
For this we have to choose xo-> which is a solution of the homogeneous Airy equa- 
tion such that it does not vanish in D a . Since the Wronskian of the couple Ai(£), 
Bi(£) is a nonzero constant, any combination with real nonzero constants of the 
form CiAi(£)+iC2Bi(£) is an everywhere nonzero solution (and the derivative is also 
nonzero). We can choose two linear independent solutions in this way, say the ones for 
which xo =Ai(£)±iBi(£). That they are formally linearly independent with respect 
to the solutions gotten from l|4.9|) follows easily, for instance from the fact that they 
correspond to different roots of the characteristic polynomial. 

To show that there is an actual solution for each formal solution in this region we 
first single out a true solution corresponding to the dominant characteristic root and 
then use it to reduce the problem to a regular one. Then we show that they give the 
expected asymptotic behavior for the solutions of the original equation. 

The ideas are similar to those used in the regular case with the exception that 
extra care is needed along the degenerate directions. 

Choose 

Xo(£)=Ai(0+ii?i(0 

and consider the nonzero formal solution that, has the leading order xo- We proceed as 
in Step 1, Section 3 to construct a rescaled recurrence with respect to the truncation 
of our formal solution. Exactly the same argument as there shows that the new 
coefficients a have smooth asymptotic series. 

What is new here is that we must provide for the estimate of the type c) to which 
end we examine the complete characteristic equation P(£, <5; A) = 0. Pisa polynomial 
in A (actually it is, to leading order, a polynomial with constant coefficients) C°° in £ 
and S. P(0,0,A) has a double root A = 1 but by assumption the second derivative does 
not vanish so that we can obtain the roots of the polynomial perturbatively. After 
series expansion, we obtain: 

i 

{ [oj(0,0) (1 + SCj + 5 2 E 3 + 0{6 3 )) + ^D x a 3 5 2 + 0(5 3 )] 

3=0 

(4.15) (1 + j X i(0S + (JX2(05 2 + j(j - l)/2x?) + 0(S 3 )) } = 

where we have taken A - 1 + Xi(O s + X2{t)$ 2 + 0(5 3 ) and 

Cj = -rOYo+Xi) 



E i = — (iVoXo + ix'iXo - JXiXo - X?) 
Xo 

Using the relations l|4.7|l we get two solutions for xi- Xi = — actually, as expected, 
we get a root Ai = 1 + 0(S S ) — and xi = 1 — ^XaXo ^ + 0{8 2 ). We see that in the 
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first order in 6 there is no root — crossing, which is not a surprize since a generic 
perturbation tends separate coalescing roots. The asymptotic series are uniformly 
valid in our domain D a . Now we show that 



(4.16) 



A 2 ((fc + l)e)-A 2 (fce) 



Xi(ke) - \ 2 (ke) 



< 



const 
lfcl + 1 



(this explains the condition c) at the begining of this section). 
We have 

1) 

| (log Xo(0)'\ > const VC + const' 
This is obvious since it holds asymptotically and the function does not vanish. Hence 

X 2 (Sk)) - Xi(Sk) > const y/5\k\ + 1 

2) 



| (log /io < const/ y/l + 1 + const' 
Using the asymptotic series for the Ai^ and the estimates 1.) and 2.) we see that 



Aa(*(fe + 1))-A a (*fc) 



X 2 (Sk) - Xi(5k) 



< 



( V£ + 1 + const' ) S 2 const 
(constVf + const') J |fc| + 1 



The similar estimates for the other roots are better but this of course does not 
improve the overall rate of convergence. Ai =1+0(S S ) and can be obviously made 
less than 1 + const S 2 in D a so that also |Ai((fe+ 1)5) — Xi(kS)\ < 1 + const' 6 2 which 
is enough for our purposes. For m ^ 1,2 we can for instance use the fact that the 
derivative of the polynomial at these points does not vanish and settle for a crude 
bound \X m (S (k + 1)) — X m (Sk)\ < const S 2 which can be obtained immediately from 
(ll~T51) . 

Now, to see that there is a true solution of the recurrence which is asymptotic to 
our formal series starting with Ai+iBi, we only have to repeat the same arguments 
as in the regular case. 



The next step is to use this particular solution to lower the order of the recurrence. 
We mimic the construction done in the proof of Lemma 13.21 to get a lower order 
recurrence in the variable du 



i-i 

Y / bs(k5,5)d k+s =0 

s=0 

(see (j^2ni) and want to check that this new recurrence satisfies the hypothesis of of 
Lemma l4.ll To leading order, the characteristic polynomial of the above equation has 
no double roots (it now has only one root equal to 1). 
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As in the regular case, &o(£>0) = — a(£, 0) = — a (0,0); fy_i(£,0) = — a;(£,0) = 
-a,(0,0). 



Noting that the coefficients of the recurrence i|3.29|l have asymptotic series valid 
throughout the domain (as finite sums of terms of the series of aj(£ tl S)Yi^+j /Yi,k), 
the boundedness of the coefficients is also trivial. 

As in the proof of Lemma 13.21 the polynomial in bj has the same roots as the poly- 
nomial in hj (except for the eliminated one) for which we have already obtained the 
estimates of type c). 

Now we have a regular problem for which we know that to each formal solution 
there is a genuine solution asymptotic to it. 



It remains to check that we can recover the asymptotic behavior of the solutions of 
the original recurrence from those of the reduced one. For the solutions corresponding 
to the characteristic roots that are less than one, exactly the same proof as in Step 3 of 
Lemma 13.21 works . For any formal solution of our original equation that corresponds 
to the largest eigenvalue and which does not vanish, the proof is the one given in 
Lemma 14.11 

In the crossing region however there might be a special interest in finding a 
particular solution which is not of exponential type and which is small for large £ 
( the Airy-like solution). To this end, a slightly different argument is necessary. We 
can obtain a formal solution of the reduced equation which is, to leading order, 



(417) Bi(£ + a)Ai(0-Bi(OAi(£ + *) 



(Ai(0-»Bi(0)(AiK + d)-*Bi(e + *)) 



(suggested by computing l|3.31|l for two formal solutions of the original equation, 
corresponding to AiztiBi)). The asymptotic behavior of (|4.17() for large £ is 



( 4 - 18 ) TIuF\ * (1 + power series) (1 « |f| « ^ 1/2 ) 

(Ai(£) - «Bi(£)r 

Writing the asymptotic representation of AiztiBi in the form 

~ C 1/A exp(2/3£ 3/2 )(l + series) 

is sufficient to see this. It is important to note that there are no powers of £ multiplying 
the asymptotics (|4.18|) ; its leading order does not vanish and (|4.18l) can be written as 
an exponential of a formal series, the form required by our arguments. 

We now apply the construction in Step 3, Section 3 to recover the solutions of 
the initial recurrence. 

Using the asymptotic behavior of the Airy functions for large argument, we get 
for the reconstructed solution the representation 



2 S ~" 4 

Y ~ £ ( ~ 1/4) exp(-£ 3 / 2 ) J2( kS ) 1/2 exp(--(fc<5)i)(l + power semes) 

3 k =J 3 
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for which the Euler-Maclaurin summation formula gives the asymptotic representa- 
tion, 



r 1/4 ex P (-^ 3/2 

for large positive £ 

In conclusion there is a true solution of the recurrence which behaves like the 
Airy function for positive k (also for negative k when this is properly interpreted) and 
our argument shows what initial conditions have to be chosen in order to obtain it. 
For negative k we see that in fact all the solutions corresponding to the largest two 
eigenvalues are comparable. 

B. The exterior region. 

Now we want to show that the solutions coming from the exterior region remain 
asymptotic to the true solutions as long as \x\ 3> e 2 ^ 3 

The problem that arises here is that the characteristic polynomial has a virtually 
degenerate root for small x and this leads to a lesser smoothness of the asymptotic 
series and ultimately to its collapse at |x| ~ e 2 / 3 . Let (5 be as in Proposition 2.2 and 
define the exterior region by 

= {x : \x\ > e' 3 } 
In what follows we make the following conventions. We write 

<?(e)=0(e°°) 

if lim c _>o e~ fc g(e) — for all k; also F{x~<) will denote a generic function such that it 
together with its derivatives of arbitrary order s satisfy the estimates 

(4.19) T{x"i)V>{x) <A S + B s x^~ s 

with the convention that |.F(a; )| = |.F(ln(a;))| < A + B\ In |sg||. 

The first step is to study the asymptotic properties of the formal solutions i.e of 
the (possibly divergent) expressions for which 

/ / oo \ 

(4.20) 2aj-(s;e)eatp(£$t(a: + ie)e*J =0 

3=0 \t=0 J 

We show by induction that <I>t and their derivatives behave like cc^t 1-4 ) and its 
derivatives. We place ourselves in the assumption of gcnericity of the crossing which 
means in particular that 9P ^' 1 ^ > const y // \x\. 

Lemma 4.3. i) If in 14-2(J\ ) the asymptotic series for the dj are of the form 
-Fkix 1 ~^ k )e k then, in the formal solution \4-2U\j we have ^t=T[xi^ 1 ^). 

ii) The same conclusion is true if aP ^' J ^ > const > uniformly in x and aj ~ 

Note that the coefficients aj of our initial recurrence are smoother then it is 
assumed in i) but this smoothness does not withstand a rescaling as done for (|3.28() . 
The proof of the Lemma is by induction on t. 
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i) It is easy to see from the eikonal equation that <&q=!F(x 3 / 2 ). Assume that the 
conclusions of the lemma hold for all t' < t. After a formal series expansion of the 
exponent of (|4.2U|) one gets 



E a jO; e ) exp 
3=0 




M s) +E 



E^ 1 E ^ s) + °(c t+1 ) 



k=0 



s > 1 

t' + S < t + 1 



(4.21) 

which, using the induction hypothesis, can be rewritten as 



E (E^V^V) exp + e* (i^+^-l*)) + 



(4.22) 



t-i 

E" A 

fc=i 



c*^ fc (a;-r- ) + 0(e t+i ) =0 



After expanding in powers of e and collecting the term in e we obtain the equation 
for $ t in the form: 



(4.23) 



(4.24) 



E 

3=0 



= 



^(z) = 



dP(X;x) 



where the derivative of the polynomial is evaluated at e = 0, A — exp(<I>Q) thus proving 
(i)- 

For (ii) the same proof works, replacing everywhere T(x x 2 ) with T{x^~ 2 ). 
Bl. Rescaling first the recurrence with respect to the approximate solution we show 
that there is a genuine solution corresponding to the maximal eigenvalue. Let ^ m be 
any function such that ty m (x; e) ~ X)o° ®m;t( x ) e and take 



(4.25) 



Y m := exp(e 1 & m (x; e)) 



The existence of a solution corresponding to the asymptotics exp(e 1 ^i(x; e)) is 
again equivalent to a solution which is 1 to all orders in e of the recurrence 



(4.26) 



^a JtS (ke,e)C k +j = 
3=0 



where ctj s (ke, e) — aj(ke,e)Y 1 ((k + j)e)/Y 1 (fce). The formal solutions of the equation 
(14.26(1 are Y m jY\. We need the roots of the new characteristic polynomial 



(4.27) 



P(X) : =E%, S (^e)A i = 
3=0 

23 



It is easy to see that the polynomial (|4.27|) has a root which is 1 to all orders in e. 
Let now G(x; e) be one of the differences G m (x; e) = v E' m (x; e) — ^i(x; e). We have, 

l 

aj, s (ke, e) exp(e^ 1 G(a; + je; e)) = o(e*) 

3=0 

for all t and so, after series expansion 
i 

(4.28) ^^(fce, £ )A J (l + ei/,(x;e)) =0(6°°) 
j=o 

where Hj(x;e) are some smooth functions of x,e and Ao := exp(Ga;(a;; e)). Using 
Lemma f4. 31 and the genericity assumptions it is not difficult to see that 

(4.29) \Hj{x;e)\ < const \x\- 1/2 

If we look for solutions of the characteristic polynomial l|4.27|l in the form Ao + 7 
we get 

1 

Q = P(X) = J2P {j) (W 
3=0 

(where the derivatives are taken with respect to A). Using (|4.28|) we obtain 7 as the 
unique small solution of the equation 



(4.30) 7 = 



which is a contraction for small enough e (and small 7 ) in the region |x| > e 13 as 
it is easy to check. We then obtain from l|4.29|l and l|4.30|) 



(4.31) I7I < const— r 

\x\ 

again valid for |x| > e 13 . 

We shall also need estimates for j(x + e) — ~f(x). Differentiating l|4.30|l with respect 
to x and using Lemma 14.31 



(4.32) \j(x + e) - j(x)\ < const e \x\~ 2 

Now we proceed as in the regular case in rewriting the recurrence in matrix form 
and evaluating the terms in the product (|3.22|l . In the matrices T := A' k 1 Ak~i the off- 
diagonal elements are estimated by T mn < const e|x| _1 and for the diagonal elements 
we have T mm = 1 + 0(e/x). Indeed, 
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X m ((p - l)e, e) - A„(f>e, e) 



A m (pe,e) - A„(pe, e) 

Since by the assumption of genericity the roots of the polynomial are separated by at 
least const y(|x|) each term in the product above can be estimated by 



1 + const e 



(A' (z) + Y(x))\x\-^ + O (e(\> (x) + Y(x))\x\- 1 / 2 ) 



< 



1 + const (ejxp 1 + e 2 \x\~ 5/2 ^ < 1 + const e\x\~ x 



in our region Ep. The nondiagonal terms are estimated in a very similar way. 

We derive the estimate ||T(a;) — I|| < Ke/x for some constant K. Assume for 
definiteness that we are on the left of the crossing point. We get, 



n r ( fc€ 



k;-ke>eP 



< const e 



-K/3 



Finally we have to control the product of the norms of the diagonal matrices 
Dfe. Since they all have one eigenvalue equal to 1 to all orders in e and for i > 2 
|Aj(x; e)J < 1 — const the only nontrivial contribution comes from A2 and this only if 
Ai and A2 have the same modulus to leading order in e. Referring to the decomposition 
A = Aq + 7 we have in this case, using Euler-Maclaurin summation formula, 

JJ \\ (ke)\=exp(e- 1 Y,m*2Uke;e)-(*iUke;e))) ~ 

fc;-fee>e' 3 \ h ) 



exp{5R($ 2;1 (- e 2 / 3 )-<i> 1;1 (- e 2 / 3 ))}< 



Also, 



n 

fe;-fee>e^ 



like) 



Ao(fce) 



< n 

fc;-fce>e' 3 



const 



\k\ 



< e 



so that also A is less than e ~ const . At this point the same arguments as in the 
regular case show that there is a true solution behaving asymptotically as the formal 
solution corresponding to the largest eigenvalue. 
B2 

Let Y\ be a solution of the recurrence relation such that Y\ ~ exp {eT l ^i;te*) 
in Ep. We now follow the same steps that led to equations l|3.28|l and (|3.29() . 

It is a matter of straightforward induction to derive from Lemma 14.31 that the 
coefficients a,- have the behavior 



k>l 
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and then clearly 

(4.33) 6,(z;e)~]>> J;fe (^-tV 

fe>0 

It is also easy to check that the recurrence (j3.2H|) is now nondegenerate in the sense 
that: 

(4.34) m[{\bo(ke)\,\b,-i(ke)\, j^j} > const > 

and the characteristic polynomial of the new recurrence does not have coalescing roots 
(the root A = 1 of (|3.28|) has been eliminated in the reduction): 

(4.35) inf{| \X m {x)\ - \X n (x)\ |} > const > (for m ^ n) 

E/3 

We are now left with a problem of the following type. Taking a recurrence of the 
form 



i 

(4.36) ^2aj(x;e)y k+ j = 

3=0 

under the following conditions: 

(4.3?) ^(^j^^^-fy 

fe>0 



(4.38) inf (|a (fce)|, |ai-i(ftc)|, . j. ,. } > const > 



(4.39) inf{| \\ m (x)\ - \X n (x)\ \ > const > (m ^ n)} 

where X m (x) are the roots of the polynomial 



i 

(4.40) P(A) := a jAke)X 3 = 

j=o 

we want to show that 

Lemma 4.4. Given a formal solution to \4.3b\j : 



exp e 1 ^Y$ t (a;) 



where §t(x) = there is a solution asymptotic to it for < x G Ep (and 

correspondingly one when x negative). 
Proof: induction on I. 
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a) We show that we can find a solution corresponding to the root the has the 
largest modulus (this will simultaneously prove the lemma for I = 1). All the ar- 
guments in Bl above apply here. Actually, now we could get some better estimates 
since we do not have small denominators in l|4.31[l , i|4.32[l and in the estimates of the 
matrices but this would not affect the final result. 

b) We assume that the conclusion is true for all recurrences of order less than 
/ — 1 and show it holds for recurrences of order By the arguments above, there is a 
true solution asymptotic to the formal solution defined by the maximum eigenvalue. 
Using it to reduce the order of the recurrence we obtain an order-Z — 1 scheme, which 
satisfies the conditions of Lemma FOI as it is easy to check and for which we thus know 
the asymptotic behavior of the solutions. It remains to verify that they can be used to 
produce solutions of the higher-order recurrence with the stated asymptotic behavior. 
For definiteness we study the subregion x < —e@. All the arguments in Step 3, Section 
3 apply if we take k\ to be — e° . The only change is that in l|3.38|) $ m; i(a;) are not 
uniformly bounded. Instead, using Lemma f4. 31 we get = .F(ln(|a;|)) < K | In e\ 
for some fixed constant K so the RHS of l (3~3gj l changes to Y m . k (l + 0{e s - 2 - K )) 

The conclusion is that after the first reduction we end up with a recurrence that 
is nondegenerate in the sense of Lemma 14.41 and for which we can control the small- 
e behavior of the solutions. Now, the reconstruction of the solutions of the original 
recurrence from the solutions of the reduced one amounts to merely repeating without 
any significant change the construction and estimates in part b) above. At this point 
in the proof it is clear that if the crossing roots are not the largest, on can reduce the 
order of the recurrence to the actual level at which the roots cross, and then apply 
the arguments above. 
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